Capítulo 5 


PROGRAMACIÓN DINÁMICA 


5.1 INTRODUCCIÓN 


Existe una serie de problemas cuyas soluciones pueden ser expresadas 
recursivamente en términos matemáticos, y posiblemente la manera más natural de 
resolverlos es mediante un algoritmo recursivo. Sin embargo, el tiempo de 
ejecución de la solución recursiva, normalmente de orden exponencial y por tanto 
impracticable, puede mejorarse substancialmente mediante la Programación 
Dinámica. 

En el diseño Divide y Vencerás del capítulo 3 veíamos cómo para resolver un 
problema lo dividíamos en subproblemas independientes, los cuales se resolvían de 
manera recursiva para combinar finalmente las soluciones y así resolver el 
problema original. El inconveniente se presenta cuando los subproblemas 
obtenidos no son independientes sino que existe solapamiento entre ellos; entonces 
es cuando una solución recursiva no resulta eficiente por la repetición de cálculos 
que conlleva. En estos casos es cuando la Programación Dinámica nos puede 
ofrecer una solución aceptable. La eficiencia de esta técnica consiste en resolver 
los subproblemas una sola vez, guardando sus soluciones en una tabla para su 
futura utilización. 


La Programación Dinámica no sólo tiene sentido aplicarla por razones de 
eficiencia, sino porque además presenta un método capaz de resolver de manera 
eficiente problemas cuya solución ha sido abordada por otras técnicas y ha 
fracasado. 


Donde tiene mayor aplicación la Programación Dinámica es en la resolución de 
problemas de optimización. En este tipo de problemas se pueden presentar distintas 
soluciones, cada una con un valor, y lo que se desea es encontrar la solución de 
valor óptimo (máximo o mínimo). 

La solución de problemas mediante esta técnica se basa en el llamado principio 
de óptimo enunciado por Bellman en 1957 y que dice: 


“En una secuencia de decisiones óptima toda subsecuencia ha de ser también 
óptima”. 


Hemos de observar que aunque este principio parece evidente no siempre es 
aplicable y por tanto es necesario verificar que se cumple para el problema en 
cuestión. Un ejemplo claro para el que no se verifica este principio aparece al tratar 
de encontrar el camino de coste máximo entre dos vértices de un grafo ponderado. 
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Para que un problema pueda ser abordado por esta técnica ha de cumplir dos 
condiciones: 


e La solución al problema ha de ser alcanzada a través de una secuencia de 
decisiones, una en cada etapa. 


e Dicha secuencia de decisiones ha de cumplir el principio de óptimo. 


En grandes líneas, el diseño de un algoritmo de Programación Dinámica consta 
de los siguientes pasos: 


1. Planteamiento de la solución como una sucesión de decisiones y verificación de 
que ésta cumple el principio de óptimo. 


2. Definición recursiva de la solución. 


3. Cálculo del valor de la solución óptima mediante una tabla en donde se 
almacenan soluciones a problemas parciales para reutilizar los cálculos. 


. Construcción de la solución óptima haciendo uso de la información contenida 
en la tabla anterior. 


ÉS 


5.2 CÁLCULO DE LOS NÚMEROS DE FIBONACCI 


Antes de abordar problemas más complejos veamos un primer ejemplo en el cual 
va a quedar reflejada toda esta problemática. Se trata del cálculo de los términos de 
la sucesión de números de Fibonacci. Dicha sucesión podemos expresarla 
recursivamente en términos matemáticos de la siguiente manera: 


E 1 sin=0,1 
AS in DA FI2) sins 


Por tanto, la forma más natural de calcular los términos de esa sucesión es 
mediante un programa recursivo: 


PROCEDURE FibRec(n: CARDINAL) : CARDINAL; 


BEGIN 
IF n<=1 THEN RETURN 1 
ELSE 
RETURN FibRec(n-1) + FibRec (n-2) 
END 
END FibRec; 


El inconveniente es que el algoritmo resultante es poco eficiente ya que su 
tiempo de ejecución es de orden exponencial, como se vió en el primer capítulo. 


Como podemos observar, la falta de eficiencia del algoritmo se debe a que se 
producen llamadas recursivas repetidas para calcular valores de la sucesión, que 
habiéndose calculado previamente, no se conserva el resultado y por tanto es 
necesario volver a calcular cada vez (véase el apartado 3.11 del capítulo 3, en 
donde se determina el número exacto de veces que se repite cada cálculo). 
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Para este problema es posible diseñar un algoritmo que en tiempo lineal lo 
resuelva mediante la construcción de una tabla que permita ir almacenando los 
cálculos realizados hasta el momento para poder reutilizarlos: 


Fib(0) | Fib(1) | Fib(2) | .. | Fib(n) 


El algoritmo iterativo que calcula la sucesión de Fibonacci utilizando tal tabla es: 


TYPE TABLA = ARRAY [O..n] OF CARDINAL 
PROCEDURE Fiblter(VAR T:TABLA;n:CARDINAL) : CARDINAL); 
VAR i:CARDINAL; 


BEGIN 
IF n<=1 THEN RETURN 1 
ELSE 
T[0] :=1; 
T[1] :=1: 


FOR i:=2 TO n DO 
T[i] :=T[i-1]+T[i-2] 
END; 
RETURN T[n] 
END 
END Fiblter; 


Existe aún otra mejora a este algoritmo, que aparece al fijarnos que únicamente 
son necesarios los dos últimos valores calculados para determinar cada término, lo 
que permite eliminar la tabla entera y quedarnos solamente con dos variables para 
almacenar los dos últimos términos: 


PROCEDURE Fiblter2(n: CARDINAL) : CARDINAL; 
VAR i,suma,x,y:CARDINAL; (* x e y son los 2 ultimos terminos *) 
BEGIN 
IF n<=1 THEN RETURN 1 
ELSE 
x:=1; y:=1; 
FOR i:=2 TO n DO 
Suma: =X+y5 Y:=X; X:=suma; 
END; 
RETURN suma 
END 
END Fiblter2; 


Aunque esta función sea de la misma complejidad temporal que la anterior 
(lineal), consigue una complejidad espacial menor, pues de ser de orden O(n) pasa 
a ser O(1) ya que hemos eliminado la tabla. 

El uso de estructuras (vectores o tablas) para eliminar la repetición de los 
cálculos, pieza clave de los algoritmos de Programación Dinámica, hace que en 
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este capítulo nos fijemos no sólo en la complejidad temporal de los algoritmos 
estudiados, sino también en su complejidad espacial. 


En general, los algoritmos obtenidos mediante la aplicación de esta técnica 
consiguen tener complejidades (espacio y tiempo) bastante razonables, pero 
debemos evitar que el tratar de obtener una complejidad temporal de orden 
polinómico conduzca a una complejidad espacial demasiado elevada, como 
veremos en alguno de los ejemplos de este capítulo. 


5.3 CÁLCULO DE LOS COEFICIENTES BINOMIALES 


En la resolución de un problema, una vez encontrada la expresión recursiva que 
define su solución, muchas veces la dificultad estriba en la creación del vector o la 
tabla que ha de conservar los resultados parciales. Así en este segundo ejemplo, 
aunque también sencillo, observamos que vamos a necesitar una tabla 
bidimensional algo más compleja. Se trata del cálculo de los coeficientes 
binomiales, definidos como: 


E) soci (=( 


El algoritmo recursivo que los calcula resulta ser de complejidad exponencial 
por la repetición de los cálculos que realiza. No obstante, es posible diseñar un 
algoritmo con un tiempo de ejecución de orden O(nk) basado en la idea del 
Triángulo de Pascal. Para ello es necesario la creación de una tabla bidimensional 
en la que ir almacenando los valores intermedios que se utilizan posteriormente: 


0 1 2 3 A k-1 k 
0 1 
1 1 1 
2 1 2 1 
3 1 3 3 1 
sl Cie) + C(nu-1,0 
SS V 


Iremos construyendo esta tabla por filas de arriba hacia abajo y de izquierda a 
derecha mediante el siguiente algoritmo de complejidad polinómica: 


PROCEDURE CoefIter(n,k: CARDINAL): CARDINAL; 
VAR 1,3]: CARDINAL; 
C: TABLA; 
BEGIN 
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FOR i:=0 TO n DO C[i,0]:=1 END; 
FOR i:=1 TO n DO C[i,1]:=i END; 
FOR i:=2 TO k DO C[i,i]:=1 END; 
FOR i:=3 TO n DO 

FOR j:=2 TO i-1 DO 


IF j<=k THEN 
C[i,j]:=C[i-1,3j-11+C[i-1,3] 
END 
END 


END; 
RETURN C[n,k] 
END CoeflIter. 


5.4 LA SUBSECUENCIA COMÚN MÁXIMA 


Hay muchos problemas para los cuales no sólo deseamos encontrar el valor de la 
solución óptima sino que además deseamos conocer cuál es la composición de esta 
solución, es decir, los elementos que forman parte de ella. En estos casos es 
necesario ir conservando no sólo los valores de las soluciones parciales, sino 
también cómo se llega a ellas. Esta información adicional puede ser almacenada en 
la misma tabla que las soluciones parciales, o bien en otra tabla al efecto. 


Veamos un ejemplo en el que se crea una tabla y a partir de ella se reconstruye 
la solución. Se trata del cálculo de la subsecuencia común máxima. Vamos en 
primer lugar a definir el problema. 


Dada una secuencia X=(x, X2 ... Xm), decimos que Z=(Z¡ Z> ... z¿) es una 
subsecuencia de X (siendo k <m) si existe una secuencia creciente (1, i> ... 14) de 
índices de X tales que para todo ¡= 1, 2, ..., k tenemos xj; = Z;. 

Dadas dos secuencias X e Y, decimos que Z es una subsecuencia común de X e Y 
si es subsecuencia de X y subsecuencia de Y. Deseamos determinar la subsecuencia 
de longitud máxima común a dos secuencias. 


Solución (O) 
Llamaremos L(i,j) a la longitud de la secuencia común máxima (SCM) de las 
secuencias X; e Y; , siendo X; el ¡-ésimo prefijo de X (esto es, X= ([x1 x> ... xp) e Y; 
el j-ésimo prefijo de Y, (Y,= (y1 y» ... y;)). 

Aplicando el principio de óptimo podemos plantear la solución como una 
sucesión de decisiones en las que en cada paso determinaremos si un carácter 
forma o no parte de la SCM. Escogiendo una estrategia hacia atrás, es decir, 


comenzando por los últimos caracteres de las dos secuencias X e Y, la solución 
viene dada por la siguiente relación en recurrencia: 


0 sii=00j=0 
Li, j)=3LG 1, ¡-D)+1 Sii0, ¡40 y x= y, 
MaxiL(i,j—1),L(i=1, JJ] sii40, ¡40 y x,4p, 
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La solución recursiva resulta ser de orden exponencial, y por tanto 
Programación Dinámica va a construir una tabla con los valores L(i, j) para evitar 
la repetición de cálculos. Para ilustrar la construcción de la tabla supondremos que 
X e Y son las secuencias de valores: 


X=(10010101)3 
Y=(010110110) 


La tabla que permite calcular la subsecuencia común máxima es: 


4 6 
1 1 


RR 00 


0 1 2 
0 


ou 
O 


3 
0 


Diag | Diag | Izq Diag 
1 Su Sup | Diag 
30 En 3 2 
Diag | Diag 
| 
] Su Sup 
] EFE 
] S Sup 
] PAE 
Diag | Diag 
| 
1 Su Sup 
] MEF 
1 Sup | Sup 
FFAA 
Diag | Diag | Su Diag | Su Diag | Sup 


Esta tabla se va construyendo por filas y rellenando de izquierda a derecha. 
Como podemos ver en cada £[i,j] hay dos datos: uno el que corresponde a la 
longitud de cada subsecuencia, y otro necesario para la construcción de la 
subsecuencia óptima. 


La solución a la subsecuencia común máxima de las secuencias X e Y se 
encuentra en el extremo inferior derecho (£[9,8]) y por tanto su longitud es seis. Si 
queremos obtener cuál es esa subsecuencia hemos de recorrer la tabla (zona 
sombreada) a partir de esta posición siguiendo la información que nos indica cómo 
obtener las longitudes óptimas a partir de su procedencia (izquierda, diagonal o 


ne 
La 


U Y 
MN ES NS) a 
e Q 
Y Y Y 
far] S Wlyg Y 
Q Q e 
a a ¿ IE 
em] 
DEE 
hb a 


¡93 


SS a) P 
pS 
al uu 
pa 


¡0/9 
[a] 
(0/9) 
Ss 


32 
U U 07) 
SN O ES 
Q 
all 


PROGRAMACIÓN DINÁMICA 183 


superior). El algoritmo para construir la tabla tiene una complejidad de orden 
O(nm), siendo n y m las longitudes de las secuencias X'e Y, 


CONSTN = ...; (* longitud maxima de una secuencia *) 

TYPE SECUENCIA = ARRAY [1..N] OF CARDINAL; 
PARES = RECORD numero: CARDINAL; procedencia:CHAR; END; 
TABLA = ARRAY [O..N],[O..N] OF PARES; 


PROCEDURE SubSecMaxima(VAR X,Y:SECUENCIA;n,m:CARDINAL;VAR L:TABLA); 
VAR i,j:CARDINAL; 
BEGIN 
FOR i:=0 TO m DO (* condiciones iniciales *) 
L[i,0] .numero:=0 
END; 
FOR j:=0 TO n DO 
L[O,j].numero:=0 
END; 
FOR i:=1 TO m DO 
FOR j:=1 TO n DO 

IF Y[il = X[j] THEN 
L[i,jl.numero:=L[i-1,j-1].numero+1; 
L[i,jl.procedencia:="D" 

ELSIF L[i-1,j].numero >=L[i,j-1].numero THEN 
L[i,jl].numero:=L[i-1,3].numero; 
L[i,jl.procedencia:="S" 

ELSE 
L[i,j].numero:=L[i,j-1].numero; 
L[i,jl.procedencia:="1I" 

END 

END 
END 
END SubSecMaxima. 


Para encontrar cuál es esa subsecuencia óptima hacemos uso de la información 
contenida en el campo procedencia de la tabla £, sabiendo que *P(por “Izq” 
significa que la información la toma de la casilla de la izquierda, *S” (“Sup”) de la 
casilla superior y de la misma manera *“D*(“Diag”) corresponde a la casilla que está 
en la diagonal. El algoritmo que recorre la tabla construyendo la solución a partir 
de esta información y de la secuencia Y es el que sigue: 


PROCEDURE Escribir(VAR L:TABLA; VAR Y:SECUENCIA; i,j:CARDINAL; 
VAR sol:SECUENCIA; VAR 1:CARDINAL); 
(* sol es la secuencia solucion, 1 su longitud, i es la longitud 
de la secuencia Y, y j la de X *) 
BEGIN 
IF (i=0) OR (¿=0) THEN RETURN END; 
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IF L[i,jl.procedencia = "D" THEN 
Escribir(L,Y,i-1,j-1,s01,1); 
sol [1] :=Y[i]; 
INC(1); 

ELSIF L[i,j].procedencia = "S" THEN 
Escribir(L,Y,i-1,j,s01,1) 

ELSE 
Escribir(L,Y,i,j-1,s01,1) 

END 

END Escribir 


La complejidad de este algoritmo es de orden O(n+m) ya que en cada paso de la 
recursión puede ir disminuyendo o bien el parámetro ¡ o bien j hasta alcanzar la 
posición £[i,j] parai=06¡=0. 


5.5 INTERESES BANCARIOS 


Dadas n funciones fi, f, ..., f, y un entero positivo M, deseamos maximizar la 
función f0c1) + f(x) +... + f(x.) sujeta a la restricción x, +x2 +... + x= M, donde 
f(0) = 0 (i=1,..,1), x, son números naturales, y todas las funciones son monótonas 
crecientes, es decir, x > y implica que f(x) > f:07). Supóngase que los valores de 
cada función se almacenan en un vector. 

Este problema tiene una aplicación real muy interesante, en donde f; representa 
la función de interés que proporciona el banco i, y lo que deseamos es maximizar el 
interés total al invertir una cantidad determinada de dinero M. Los valores x, van a 
representar la cantidad a invertir en cada uno de los n bancos. 


Solución (O) 


Sea f; un vector que almacena el interés del banco i (1 < ¡ < n) para una inversión 
de 1, 2, 3, .... M pesetas. Esto es, f(7) indicará el interés que ofrece el banco ¡ para ¡ 
pesetas, con0<i¡<n,0<j<M. 

Para poder plantear el problema como una sucesión de decisiones, llamaremos 
[,(M) al interés máximo al invertir M pesetas en n bancos, 


L(M) = 4100) +2) + ... + ful) 
que es la función a maximizar, sujeta a la restricción x¡ +x2+... +x,=M. 
Veamos cómo aplicar el principio de óptimo. Si /,(M) es el resultado de una 
secuencia de decisiones y resulta ser óptima para el problema de invertir una 


cantidad M en n bancos, cualquiera de sus subsecuencias de decisiones ha de ser 
también óptima y así la cantidad 


Ii M— xp) =f1001) + fa) +... + fa 101) 


será también óptima para el subproblema de invertir (M — x,) pesetas en n — 1 
bancos. Y por tanto el principio de óptimo nos lleva a plantear la siguiente relación 
en recurrencia: 
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AE) sin=1 


[, (x) = Maxtl,_(x=0)+ f,(0)). en otro caso. 
<t<x 


[5.1] 


Para resolverla y calcular /,(M), vamos a utilizar una matriz / de dimensión nxM 
en donde iremos almacenando los resultados parciales y así eliminar la repetición 
de los cálculos. El valor de /[1,/] va a representar el interés de j pesetas cuando se 
dispone de ¡ bancos, por tanto la solución buscada se encontrará en /[n,M]. Para 
guardar los datos iniciales del problema vamos a utilizar otra matriz F, de la misma 
dimensión, y donde F[i,¡/] representa el interés del banco ¡ para j pesetas. 


En consecuencia, para calcular el valor pedido de /[n,M] rellenaremos la tabla 
por filas, empezando por los valores iniciales de la ecuación en recurrencia, y 
según el siguiente algoritmo: 


CONSTn ...; (tk numero de bancos x*) 
M=...; (tk cantidad a invertir x*) 
TYPE MATRIZ = ARRAY [1..n],[O..M] OF CARDINAL; 


PROCEDURE Intereses(VAR F:MATRIZ;VAR 1I:MATRIZ) : CARDINAL; 
VAR i,j: CARDINAL; 
BEGIN 
FOR i:=1 TO n DO 1[i,0]:=0 END; 
FOR j:=1 TO M DO 1[1,3j]:=F[1,3] END; 
FOR i:=2 TO n DO 
FOR j:=1 TO M DO 
1[1,3]:=Max(1,F,i,j) 
END 
END; 
RETURN 1[n,M] 
END Intereses; 


La función Max es la que calcula el máximo que aparece en la expresión [5.1]: 


PROCEDURE Max(VAR 1,F:MATRIZ;i,3j:CARDINAL) : CARDINAL; 
VAR max,t: CARDINAL; 
BEGIN 
max:= I[i-1,j] + F[i,01; 
FOR t:=1 TO j DO 
max :=Max2 (max, I[i-1,j-t]+F[i,t]) 
END; 
RETURN max 
END Max; 


La función Max2 es la que calcula el máximo de dos números naturales. La 
complejidad del algoritmo completo es de orden O(nM?), puesto que la 
complejidad de Max es O(¡) y se invoca dentro de dos bucles anidados que se 
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desarrollan desde 1 hasta M. Es importante hacer notar el uso de parámetros por 
referencia en lugar de por valor para evitar la copia de las matrices en la pila de 
ejecución del programa. 


Por otro lado, la complejidad espacial del algoritmo es del orden O(1M), pues 
de este orden son las dos matrices que se utilizan para almacenar los resultados 
intermedios. 


En este ejemplo queda de manifiesto la efectividad del uso de estructuras en los 
algoritmos de Programación Dinámica para conseguir obtener tiempos de 
ejecución de orden polinómico, frente a los tiempos exponenciales de los 
algoritmos recursivos iniciales. 


5.6 EL VIAJE MÁS BARATO POR RÍO 


Sobre el río Guadalhorce hay n embarcaderos. En cada uno de ellos se puede 
alquilar un bote que permite ir a cualquier otro embarcadero río abajo (es imposible 
ir río arriba). Existe una tabla de tarifas que indica el coste del viaje del 
embarcadero í al j para cualquier embarcadero de partida ¡ y cualquier embarcadero 
de llegada ¡ más abajo en el río (i < 7). Puede suceder que un viaje de ¡ a j sea más 
caro que una sucesión de viajes más cortos, en cuyo caso se tomaría un primer bote 
hasta un embarcadero k y un segundo bote para continuar a partir de k. No hay 
coste adicional por cambiar de bote. 


Nuestro problema consiste en diseñar un algoritmo eficiente que determine el 
coste mínimo para cada par de puntos i,j (i < j) y determinar, en función de n, el 
tiempo empleado por el algoritmo. 


Solución (O) 


Llamaremos 7]i,j] a la tarifa para ir del embarcadero ¡ al j (directo). Estos valores 
se almacenarán en una matriz triangular superior de orden n, siendo n el número de 
embarcaderos. 


El problema puede resolverse mediante Programación Dinámica ya que para 
calcular el coste óptimo para ir del embarcadero i al ¡ podemos hacerlo de forma 
recurrente, suponiendo que la primera parada la realizamos en un embarcadero 
intermedio k (1<k< J): 

C(,j) = T(1,k) + C(k,J). 


En esta ecuación se contempla el viaje directo, que corresponde al caso en el 
que k coincide con j. Esta ecuación verifica también que la solución buscada C(1,j) 
satisface el principio del óptimo, pues el coste C(k,j), que forma parte de la 
solución, ha de ser, a su vez, óptimo. Podemos plantear entonces la siguiente 
expresión de la solución: 


0 sii=j 


CODE MTM) +C (Y sis] 152] 


La idea de esta segunda expresión surge al observar que en cualquiera de los 
trayectos siempre existe un primer salto inicial óptimo. 
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Para resolverla según la técnica de Programación Dinámica, hace falta utilizar 
una estructura para almacenar resultados intermedios y evitar la repetición de los 
cálculos. La estructura que usaremos es una matriz triangular de costes Cl[i,j], que 
iremos rellenando por diagonales mediante el procedimiento que hemos 
denominado Costes. La solución al problema es la propia tabla, y sus valores Cl[i,7] 
indican el coste óptimo para ir del embarcadero ¡al 7. 


CONST MAXEMBARCADEROS = ...; 
TYPE MATRIZ=ARRAY [1. .MAXEMBARCADEROS] , [1. .MAXEMBARCADEROS] OF CARDINAL); 


PROCEDURE Costes(VAR C:MATRIZ;n:CARDINAL) ; 
VAR i, diagonal: CARDINAL; 
BEGIN 
FOR i:=1 TO n DO C[i,i]:=0 END; (* condiciones iniciales x) 
FOR diagonal:=1 TO n-1 DO 
FOR i:=1 TO n-diagonal DO 
C[i,i+diagonal]:=Min(C,i,i+diagonal) 
END 
END 
END Costes; 


Dicho procedimiento utiliza la siguiente función, que permite calcular la 
expresión del mínimo que aparece en la ecuación en recurrencia [5.2]: 


PROCEDURE Min(VAR C:MATRIZ; i,j:CARDINAL): CARDINAL; 
VAR k,min:CARDINAL; 
BEGIN 
min:=MAX (CARDINAL) ; 
FOR k:=i+1 TO j DO 
min:=Min2(min,T[i,k] + C[k,j31) 
END; 
RETURN min 
END Min; 


La función Min2 es la que calcula el mínimo de dos números naturales. Es 
importante observar que esta función, por la forma en que se va rellenando la 
matriz C, sólo hace uso de los elementos calculados hasta el momento. 


La complejidad del algoritmo es de orden O(n”), pues está compuesto por dos 
bucles anidados de tamaño n, que contienen la llamada a una función de orden 
O(n), la que calcula el mínimo. 


5.7 TRANSFORMACIÓN DE CADENAS 


Sean u y v dos cadenas de caracteres. Se desea transformar u en v con el mínimo 
número de operaciones básicas del tipo siguiente: eliminar un carácter, añadir un 
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carácter, y cambiar un carácter. Por ejemplo, podemos pasar de abbac a abcbc en 
tres pasos: 


abbac —> abac (eliminamos b en la posición 3) 
> ababc (añadimos b en la posición 4) 
> abcbc (cambiamos a en la posición 3 por c) 


Sin embargo, esta tranformación no es óptima. Lo que queremos en este caso es 
diseñar un algoritmo que calcule el número mínimo de operaciones, de esos tres 
tipos, necesarias para transformar u en v y cuáles son esas operaciones, estudiando 
su complejidad en función de las longitudes de u y v. 


Solución (8) 


En primer lugar, la transformación mostrada arriba no es óptima ya que podemos 
pasar de abbac a abcbc en sólo dos pasos: 


abbac —> abcac (cambiamos b en la posicion 3 por c) 
> abcbc (cambiamos a en la posición 4 por c) 


Llamaremos m a la longitud de la cadena u, n a la longitud de la cadena v, y 
OB(m,n) indicará el número de operaciones básicas mínimo para transformar una 
cadena u de longitud m en otra cadena v de longitud ». 

Para resolver el problema utilizando Programación Dinámica es necesario 
plantearlo como una sucesión de decisiones que satisfaga el principio de óptimo. 

Para plantearla, vamos a fijarnos en el último elemento de cada una de las 
cadenas. Si los dos son iguales, entonces tendremos que calcular el número de 
operaciones básicas necesarias para obtener de la primera cadena menos el último 
elemento, y la segunda cadena también sin el último elemento, es decir, 


OB(m,n) = OB(m-1,n—-1) Si 4, = Vp. 


Pero si los últimos elementos fueran distintos habría que escoger la situación 
más beneficiosa de entre tres posibles: (i) considerar la primera cadena y la 
segunda pero sin el último elemento, o bien (ii) la primera cadena menos el último 
elemento y la segunda cadena, o bien (iii) las dos cadenas sin el último elemento. 
Esto da lugar a la siguiente relación en recurrencia para OB(m,n) para este caso: 


OB(m,nyF1+MinfOB(m,n—1),OB(m-1,n),OB(m-1,n-1)j si mA, nA y UV. 
En cuanto a las condiciones iniciales, tenemos las tres siguientes: 

OB(0,0)=0, OB(m,0)=m y OB(0,n)=n. 
Una vez disponemos de la ecuación en recurrencia necesitamos resolverla 


utilizando alguna estructura que nos permita reutilizar resultados intermedios, 
como mostramos a continuación: 


CONST MAXCARACTERES = ...; 
TYPE CADENA=ARRAY [1. .MAXCARACTERES] OF CHAR; 
TABLA=ARRAY [O. .MAXCARACTERES] , [0. .MAXCARACTERES] OF CARDINAL; 
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PROCEDURE Cadena(VAR OB:TABLA;u,v:CADENA;n,m: CARDINAL) : CARDINAL; 
VAR 1,3: CARDINAL; 
BEGIN 
FOR i:=0 TO m DO OB[i,0]:=1; END; 
FOR j:=0 TO n DO OB[O,j]:=j END; 
FOR i:=1 TO m DO 
FOR j:=1 TO n DO 
IF u[lil= v[j] THEN 0B[i,j3]:=0B[i-1,j-1] 
ELSE 0B[i,3]:=Min3(0B[i,3j-11,0B[i-1,3],0B[i-1,j-11)+1; 
END 
END 
END; 
RETURN OB[m,n] 
END Cadena; 


El procedimiento Cadena va a permitir la creación de la tabla OB que calcula el 
número mínimo de operaciones básicas. La solución se encuentra en OB[m,n], y la 
tabla se construye fila a fila (a partir de los valores que definen las condiciones 
iniciales) para poder ir reutilizando los valores calculados previamente. La función 
Min3 es la que calcula el mínimo de tres enteros. 


Como el algoritmo se limita a dos bucles anidados que sólo incluyen 
operaciones constantes la complejidad de este algoritmo es de orden O(mn). 


5.8 LA FUNCIÓN DE ACKERMANN 


La función de Ackermamn se define recursivamente del modo siguiente: 


Ack(0,n)=n>+1 
Ack (m,0) = Ack(m—1,1) sim >0 
Ack(m,n)= Ack(m —1, Ack(m,n—1)) si m,n>0. 


Nos planteamos los beneficios de diseñar, si es posible, un algoritmo de 
Programación Dinámica para calcular Ack(m,n). 


Solución (O) 


Este problema nos permite analizar uno de los aspectos más importantes de la 
Programación Dinámica: la búsqueda de estructuras que permitan resolver una 
ecuación en recurrencia reutilizando los cálculos realizados hasta el momento. 
Como veremos en este ejemplo, esta tarea es a veces complicada. 


Si observamos la definición recursiva de esta función para valores de m y n 
mayores que 0, vemos que para calcular el valor de Ack(m,n) será necesario utilizar 
un vector 4 suficientemente grande sobre el cual se irán almacenando de izquierda 
a derecha los sucesivos valores de la función, comenzando con m = 0. En cada paso 
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m del algoritmo se actualizarán los elementos A[i] (1 < ¡ < n), que van a almacenar 
el valor de Ack(m,i). 

Para el cálculo de cada elemento A[i] se necesitará, además del elemento 
anterior (obsérvese que A|i-1] contiene el valor de Ack(m,i-1)), un elemento del 
vector calculado en el paso anterior. La dificultad que entraña es que, conforme 
aumenta m, el elemento al que tenemos que referirnos del vector previamente 
calculado tiene un índice excesivamente elevado, 


2 


con lo cual la dimensión del vector 4 ha de ser muy grande. 


Un algoritmo que resuelve el problema siguiendo estas indicaciones es el 
siguiente: 


9.2” 


CONST MaxIndice = ...; 
TYPE VECTOR = ARRAY[O..MaxIndice] OF CARDINAL; 


PROCEDURE Ackerman(m,n:CARDINAL): CARDINAL; 
VAR i,j:CARDINAL; max,1:LONGREAL; A: VECTOR; 
BEGIN 
max:=1.0; 
FOR i:=1 TO m DO 
max:=Pot(2.0,max) 


END; 

FOR i:=0 TO VAL(CARDINAL,Pot (max ,LONGREAL(n))) DO 
A[i] :=i+1 

END; 

1:=max; 

FOR i:=1 TO m DO 
A[oJ:=A[1]; 


1:=Log(2.0,1); 
FOR j:=1 TO VAL(CARDINAL,Pot(1,LONGREAL(n))) DO 
A[j]:=A[A[5-11] 
END; 
END; 
RETURN A[n] 
END Ackerman; 


La función Pot(a,b) es la que calcula la potencia b-ésima de un número a dado, 
esto es, a”, y la función Log(a,b) la que calcula el logaritmo en base a de b. 


La complejidad temporal del algoritmo viene determinada en primer lugar por 
el valor del parámetro m ya que ha de actualizarse m veces el vector 4, y además 
por el tamaño de este vector, resultando en un orden O(m:-MaxIndice) debido a los 
dos bucles anidados del programa. 

Este procedimiento es, sin embargo, muy “ingenuo”. Y decimos esto porque, 
aunque perfectamente correcto desde un punto de vista teórico, su utilidad práctica 
es más bien poca. Al tener que manejar números tan grandes, que a su vez deben 
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ser usados como índices del vector, es muy pequeño el número de pasos que 
soporta sin exceder la capacidad de cálculo de cualquier ordenador. 
Desgraciadamente no conseguimos de esta forma manejar la “intratabilidad” de la 
función de Ackerman. 


5.9 EL PROBLEMA DEL CAMBIO 


Dentro del tema dedicado a algoritmos ávidos vimos un algoritmo para minimizar, 
dado un sistema monetario, el número de monedas necesarias para reunir una 
cantidad. Aquel algoritmo funcionaba cuando los tipos de monedas eran, por 
ejemplo, de 1, 5, 10 y 25 unidades, pero no obtenía necesariamente la 
descomposición óptima si añadíamos una moneda de 12 unidades al sistema. 


Dado que el algoritmo ávido para este problema falla en algunas ocasiones, nos 
planteamos si puede resolverse utilizando Programación Dinámica de forma que la 
solución sea satisfactoria en todos los casos. 


Solución (8) 


Sea n el número de tipos de monedas distintos, L la cantidad a conseguir y 7]1..n] 
un vector con el valor de cada tipo de moneda del sistema. Supondremos que 
disponemos de una cantidad inagotable de monedas de cada tipo. 


Llamaremos C(i,j) (1 <i<nm, 1 << £L) al número mínimo de monedas para 
obtener la cantidad j restringiéndose a los tipos 7|1], 7]2], ..., 7[1]. Si no se puede 
conseguir dicha cantidad entonces C(i,j) = co, En primer lugar hemos de encontrar 
una expresión recursiva de C(1,/). Para ello observemos que en cada paso existen 
dos opciones: 


1. No incluir ninguna moneda del tipo 7(i). Esto supone que el valor de C(i,/) va a 
coincidir con el de C(¡-1,7), y por tanto C(1,j) = C(i-1,]). 

2. Sí incluirla. Pero entonces, al incluir la moneda del tipo 7(1), el número de 
monedas global coincide con el número óptimo de monedas para una cantidad 
(¡ - T(1)) más esta moneda 7(1) que se incluye, es decir podemos expresar C(1,j) 
en este caso como C(i,j)=1+C(,j-T(0). 
El cálculo de C(,7) óptimo tomará la solución más favorable, es decir, el menor 

valor de ambas opciones. Con esto, la relación en recurrencia queda definida como: 


00 sii=1y1<j¡<T(1) 

0 sij=0 
C(i,j)=31+C(1,j-T()) sii=1ly ¡2T(i) 

C(i—1, j) sii>lyj<T() 


MiniC(¡=1, 5),1+C(i, ¡-T())) en otro caso 


Una vez disponemos de la solución recursiva del problema, aplicaremos un 
algoritmo de Programación Dinámica para calcular los C(n,j), 1 <j < L, mediante 
el uso de un vector de longitud £. 
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Llamemos C a dicho vector, que verifica que en cada paso ¡(1 <i¿<n), Cl¡] vaa 
contener el valor de C(i,7). La idea es ir actualizando dicho vector paso a paso hasta 
llegar al paso n. El algoritmo que construye este vector es el siguiente: 


CONST n .; Gr num. tipos de monedas distintos del sistema *) 
L=...; (* cantidad a conseguir *) 

TYPE  TIPOMONEDA = ARRAY[1..n] OF CARDINAL; 
VECTOR = ARRAY[O..L] OF CARDINAL; 


PROCEDURE Cambio(VAR C:VECTOR;L,n:CARDINAL;VAR T:TIPOMONEDA) : CARDINAL; 
VAR 1,j:CARDINAL; 
BEGIN 
C[o] :=0; 
FOR i:=1 TO n DO 
FOR j:=1 TO L DO 
IF (i=1) AND (¿<T[i]) THEN 
C[j] :=MAX (CARDINAL) 
ELSIF i=1 THEN 
C[j]:=1+C[5-T[111] 
ELSIF j>=T[i] THEN 
C[jl:=Min2(C[j],1+C[j-T[i11) 
(* ELSE C[jl] no se modifica *) 
END 
END; 
END; 
RETURN C[L] 
END Cambio; 


El algoritmo devuelve un valor, que es el número óptimo de monedas necesario 
para obtener la cantidad L, siendo su complejidad temporal de orden O(n£L) y la 
espacial de orden O(£). 


Una vez calculado el vector que nos permite encontrar el número mínimo de 
monedas es posible diseñar un algoritmo que construya la solución. 


Para ello va a ser necesario mantener una tabla de valores lógicos P[1,/] que 
indique la procedencia del valor C(i,/) en la expresión en recurrencia, es decir, si 
C(,j) = CG—1,j) (lo que indica que no se toma una moneda de valor 7[1]) o bien 
C(,)= CG j-T[1]) + 1 (indicando que sí se incluye una moneda del valor 7]i]). Por 
tanto, bastará definir P[i,¡] = FALSE en el primer caso, TRUE en el segundo. 

Para poder calcular los valores de esta matriz de procedencia P[i,/] se necesitan 
tener presentes todos los valores de C(i,), para lo cual ya no es suficiente un vector 
C como el utilizado en el apartado anterior, sino que será necesaria la creación de 
una matriz Cli,/] que conserve los distintos valores para ¡ = 1,2,...,n. El algoritmo 
que implementa tal estrategia es el siguiente: 


TYPE MONEDAS = ARRAY[1..n] OF CARDINAL; 
MATRIZ = ARRAY[1..n],[0..L] OF BOOLEAN; 
CAMBIO = ARRAY[1..n],[0..L] OF CARDINAL; 
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PROCEDURE Procedencia(VAR P:MATRIZ; VAR C:CAMBIO; L,n:CARDINAL; 
T:MONEDAS) ; 
VAR 1,j:CARDINAL; 
BEGIN 
FOR i:=1 TO n DO 
P[i,0] :=FALSE; 
C[i,0]:=0 
END; 
FOR i:=1 TO n DO 
FOR j:=1 TO L DO 
IF (i=1) AND (¿<T[i]) THEN 
C[i,j]:=MAX(CARDINAL) ; 
P[i,jl:=FALSE 
ELSIF i=1 THEN 
Cli,jl:=1 + Cli,j-T[11]; 
P[i,j]:=TRUE 
ELSIF ¿¡<T[il THEN 
Cli,j1:=C[i-1,3]; 
P[i,jl:=FALSE 
ELSE 
C[i,j]:=Min2(C[i-1,3]1,1+C[i,j-T[i11); 
P[i,j1:=(C[i,j] <> Cli-1,3]) 
END 
END 
END 
END Procedencia; 


La solución se encuentra recorriendo la tabla P en sentido inverso, comenzando 
por el valor P[n, £], como se muestra en el siguiente algoritmo: 


TYPE NUMMONEDAS = ARRAY[1..n] OF CARDINAL; 


PROCEDURE Monedas (P:MATRIZ;C:CAMBIO;L,n:CARDINAL;T:MONEDAS) : NUMMONEDAS; 
VAR monedas : NUMMONEDAS; ind,i,j:CARDINAL; 
BEGIN 
i:=n; j:=L; 
FOR ind:=1 TO n DO monedas lind] :=0 END; 
WHILE (i<>0) AND (¿<>0) DO 
IF P[i,j]l=FALSE THEN DEC(i) 
ELSE monedas[il]:=monedas [1]+1; j:=j-T[i] 
END 
END; 
IF i=0 THEN monedas [1] :=C[i,j]+monedas[1] END; 
RETURN monedas 
END Monedas; 
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La resolución del problema requiere por una parte el algoritmo Procedencia 
para la construcción de la tabla P que permite conocer el número de monedas, y 
por tanto cuando se trata de n tipos de monedas diferentes y una cantidad L, su 
orden de complejidad es O(n£). 


Además, para conocer la solución es necesario recorrer la tabla desde la 
posición P[n,£] hasta P[0,0] a través de n — 1 pasos de orden de complejidad O(n) 
para llegar a la fila 1. Asimismo hay que tener en cuenta los pasos que hay que dar 
de derecha a izquierda hasta llegar a la columna cero, que viene dado por el 
número de monedas que intervienen en la solución, es decir, por el valor de C|n,£]. 
Podemos concluir por tanto que su complejidad es O(n+C[n,£]). 


5.10 EL ALGORITMO DE DIJKSTRA 


Sea un grafo ponderado g = (V,4), donde Y es su conjunto de vértices, A el 
conjunto de arcos y sea £Lli,j] su matriz de adyacencia. Queremos calcular el 
camino más corto entre un vértice v, tomado como origen y cada vértice restante v; 
del grafo. 


El clásico algoritmo de Dijkstra trabaja en etapas, en donde en cada una de ellas 
va añadiendo un vértice al conjunto D que representa aquellos vértices para los que 
se conoce su distancia al vértice origen. Inicialmente el conjunto D contiene sólo al 
vértice origen. 


Aún siendo el algoritmo de Dijkstra un claro ejemplo de algoritmo ávido, nos 
preguntamos si puede ser planteado como un algoritmo de Programación 
Dinámica, y si de ello se deriva alguna ventaja. 


Solución (O) 


La técnica de la Programación Dinámica tiene grandes ventajas, y una de ellas es la 
de ofrecer un diseño adecuado y eficiente a todos los problemas que puedan 
plantearse de forma recursiva y cumplan el principio del óptimo. 


Así, es posible plantear el algoritmo de Dijkstra en términos de la Programación 
Dinámica, y de esta forma aprovechar el método de diseño y las ventajas que esta 
técnica ofrece. 


En primer lugar, observemos que es posible aplicar el principio de óptimo en 
este caso: si en el camino mínimo de v, a v;, está un vértice v, como intermedio, los 
caminos parciales de v,a v, y de v, a v; han de ser a su vez mínimos. 


Llamaremos D(¡) al vector que contiene el camino mínimo desde el vértice 
origen ¡= 1 a cada vértice v;, 2 <¡< n, siendo n el número de vértices. Inicialmente 
D contiene los arcos L(1,j), o bien oo si no existe el arco. A continuación, y para 
cada vértice v, del grafo con k % 1, se repetirá: 


D(j) =MinfDJ), D(k) + L(k, DY [5.3] 


De esta forma el algoritmo que resuelve el problema puede ser implementado 
como sigue: 
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CONST n = ...; (* numero de vertices del grafo *) 
TYPE MATRIZ = ARRAY [1..n],[1..n] OF CARDINAL; 
MARCA = ARRAY [1..n] OF BOOLEAN; (* elementos ya consideradosx*) 


SOLUCION = ARRAY [2..n] OF CARDINAL; 


PROCEDURE Dijkstra(VAR L:MATRIZ;VAR D: SOLUCION); 
VAR i,j,menor,pos,s:CARDINAL; S:MARCA; 
BEGIN 
FOR i:=2 TO n DO 
S[i]:=FALSE; 
D[i]:=L [1,1] 
END; 
S[1]:=TRUE; 
FOR i:=2 TO n-1 DO 
menor :=Menor (D,S,pos) ; 
S[pos] :=TRUE; 
FOR j:=2 TO n DO 
IF NOT(S[j]1) THEN 
D[j]:= Min20D[j],D[pos]+L[pos,j3]) 
END; 
END; 
END 
END Dijkstra; 


La función Menor es la que calcula el mínimo de la expresión en recurrencia 
[5.3] que define la solución del problema: 


PROCEDURE Menor(VAR D:SOLUCION; VAR S:MARCA; VAR pos:CARDINAL) 
: CARDINAL; 
VAR menor,i:CARDINAL; 
BEGIN 
menor :=MAX(CARDINAL); pos:=1; 
FOR i:=2 TO n DO 
IF NOT(S[i]) THEN 
IF D[i]<menor THEN 
menor :=D[i]; pos:=i 
END 
END 
END; 
RETURN menor 
END Menor; 


La complejidad temporal del algoritmo es de orden O(x?), siendo de orden O(n) 
su complejidad espacial. No ganamos sustancialmente en eficiencia mediante el 
uso de esta técnica frente al planteamiento ávido del algoritmo, pero sin embargo sí 
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ganamos en sencillez del diseño e implementación de la solución a partir del 
planteamiento del problema. 


5.11 EL ALGORITMO DE FLOYD 


Sea g un grafo dirigido y ponderado. Para calcular el menor de los caminos 
mínimos entre dos vértices cualesquiera del grafo, podemos aplicar el algoritmo de 
Dijkstra a todos los pares posibles y calcular su mínimo, o bien aplicamos el 
siguiente algoritmo (Floyd) que, dada la matriz L de adyacencia del grafo g, calcula 
una matriz D con la longitud del camino mínimo que une cada par de vértices: 


CONST n = ...; (* numero de vertices del grafo *) 
TYPE MATRIZ = ARRAY[1..n],[1..n] OF CARDINAL; 


PROCEDURE Floyd (VAR L,D:MATRIZ); 
VAR 1,],k: CARDINAL; 
BEGIN 
FOR i:=1 TO n DO FOR j:=1 TO n DO 
D[i,j]:=L[i,3] 
END END; 
FOR k:=1 TO n DO 
FOR i:=1 TO n DO 
FOR j:=1 TO n DO 
D[i,j]:=Min20D[i,3],D[i,k]+D[k,3]) 
END 
END 
END 
END Floyd; 


Nos planteamos si tal algoritmo puede ser considerado o no de Programación 
Dinámica, es decir, si reune las características esenciales de ese tipo de algoritmos. 


Solución (O) 


Este algoritmo puede ser considerado de Programación Dinámica ya que es 
aplicable el principio de óptimo, que puede enunciarse para este problema de la 
siguiente forma: si en el camino mínimo de v;a v;, ves un vértice intermedio, los 
caminos de v; a v, y de v¿a v; han de ser a su vez caminos mínimos. Por lo tanto, 
puede plantearse la relación en recurrencia que resuelve el problema como: 


D,(i,j)=MintD, ,(0,Dy_, (de J)) 


Tal ecuación queda resuelta mediante el algoritmo presentado que, siguiendo el 
esquema de la Programación Dinámica, utiliza una matriz para evitar la repetición 
de los cálculos. Con ello consigue que su complejidad temporal sea de orden O(n*) 
debido al triple bucle anidado en cuyo interior hay tan sólo operaciones constantes. 
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5.12 EL ALGORITMO DE WARSHALL 


Al igual que ocurre con el algoritmo de Floyd descrito en el apartado anterior, 
estamos interesados en encontrar caminos entre cada dos vértices de un grafo. Sin 
embargo, aquí no nos importa su longitud, sino sólo su existencia. Por tanto, lo que 
deseamos es diseñar un algoritmo que permita conocer si dos vértices de un grafo 
están conectados o no, lo que nos llevaría al cierre transitivo del grafo. 


Solución (O) 


Para un grafo g = (VA) cuya matriz de adyacencia sea £, el algoritmo pedido puede 
ser implementado como sigue: 


CONST n = ...; (* numero de vertices del grafo *) 
TYPE MATRIZ = ARRAY[1..n],[1..n] OF BOOLEAN; 


PROCEDURE Warshal1l (VAR L,D:MATRIZ); 
VAR i,j,k: CARDINAL; 
BEGIN 
FOR i:=1 TO n DO 
FOR j:=1 TO n DO 
D[i,j1:=L[i,3] 
END 
END; 
FOR k:=1 TO n DO 
FOR i:=1 TO n DO 
FOR j:=1 TO n DO 
D[i,jl1:=D[i,j] OR (D[i,k] AND D[x,3]) 
END 
END 
END 
END Warshall; 


Tras la ejecución del algoritmo, la solución se encuentra en la matriz D, que 
verifica que Di] = TRUE si y sólo si existe un camino entre los vértices i y j¿. 
Obsérvese la similitud entre este algoritmo y el de Floyd. En cuanto a su 
complejidad, podemos afirmar que es de orden O(n 3) debido al triple bucle anidado 
que posee, en cuyo interior sólo se realizan operaciones constantes. 


5.13 ORDENACIONES DE OBJETOS ENTRE DOS RELACIONES 


Dados n objetos, queremos calcular el número de ordenaciones posibles según las 
relaciones “<” e “=”. Por ejemplo, dados tres objetos A, B y C, el algoritmo debe 
determinar que existen 13 ordenaciones distintas: A=B=C, A=B<C, A<B=C, 
A<B<C, A<C<B, A=C<B, B<A=C, B<A<C, B<C<A, B=C<A, C<A=B, C<A<B y 
C<B<A. 
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Solución (és) 


Llamaremos C, al número de ordenaciones posible con n objetos. Si 1 <k< n, 
podemos expresar C, como: 
k-1 
C, =19419 +. +19 =Y 10 


¡=0 
siendo / o el número de formas posibles de poner k elementos en donde hay j 


6» 


símbolos (es decir, k—¡— 1 símbolos distintos), con 0 << k, 1 <k<n. Vamos 
a tratar de expresar C; en función de C;.¡. Para ello, supongamos que ya tenemos 


k2 
== (k-1) 
Cu= 1 
¡>0 


y añadimos un nuevo elemento. Pueden ocurrir dos casos: que sea distinto a todos 
los k-1 elementos anteriores, o bien que sea igual a uno de ellos. Entonces, la 
expresión de C; va a venir dada por: 


Co = ALE + (DL + (DY 2) + 
k-2 

(e DIES + (DA 4. 4 UE ALEDO DIA 
j=0 


LE k-1 z 
Con esto, tenemos C; en función de los 14 le es decir, de los componentes del 
caso anterior. Ahora bien, es posible también relacionar los 7 Peon los 1” de la 


siguiente manera: 
LO = 1000 
7 =(k- E + (k-— NES 
Pia =(k-— 2 + (k — 2 


4) 9700 9700 
PDA DAME, 3 E 
(4) — pun 
Mrs Ss Lis 


cuyas condiciones iniciales son / E =2 4 0 =1. Esto también puede expresarse 


como sigue: 
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IP =(k= UA +I EP) para0< ¡<k-1ly2<k<n 


ds 
| 
1? =0 
e =0 


Así, el problema puede resolverse calculando cada E (0 < 7 < n—1), para 


finalmente calcular C, mediante la expresión: 


n—-1 
=> (1) 
E 
520 
El algoritmo que implementa tal estrategia es el siguiente: 


CONST n = ...; (* numero de objetos *) 
TYPE VECTOR = ARRAY [-1..n] OF INTEGER; 


PROCEDURE Ordenaciones(VAR I:VECTOR; n: INTEGER) : INTEGER; 
VAR x,y,S,k,]: INTEGER; 
BEGIN 
IF n <= 1 THEN RETURN n END; (* caso base *) 
FOR j:=-1 TO n DO 1[j]:=0 END; (* inicializamos el vector 1 *) 
I[0] :=1; x:=0; 
FOR k:=2 TO n DO 
FOR j:=0 TO n-1 DO 
IF j>1 THEN 1[j-2]:=y END; 
y:=x; 
x:=(k-3)*(1[31+1[3-11) 
END; 
I[n-2] :=y; I[n-1):=x; 
END; 
s:=0; 
FOR j:=0 TO n-1 DO 
s:=s+I[j] 
END; 
RETURN s 
END Ordenaciones; 


Respecto a su complejidad espacial, tan sólo utiliza el vector / por lo que es de 
orden O(n). Y en cuanto a su complejidad temporal, el algoritmo utiliza dos bucles 
anidados para el cálculo de los valores del vector, por lo que podemos afirmar que 
su orden es O(n?). 
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5.14 EL VIAJANTE DE COMERCIO 


¿Podría aplicarse la técnica de Programación Dinámica al problema del viajante de 
comercio? Recordemos que este problema consistía en encontrar el camino sin 
ciclos de menor peso de un grafo ponderado que recorra todos los vértices y vuelva 
al vértice original. 


Solución (O) 


En primer lugar, vamos a plantear la solución del problema como una sucesión de 
decisiones que verifique el principio de óptimo. La idea va a consistir en construir 
una solución mediante la búsqueda sucesiva de recorridos mínimos de tamaño 1, 2, 
3, etc. 


Representando el problema a través de un grafo g = (V,4), y siendo L su matriz 
de adyacencia, cada recorrido del viajante que parte del vértice v; estará formado 
por un arco (v,,v;) para algún vértice v, perteneciente a V-fv,) y un camino de vz al 
vértice v;. 

Pero si el recorrido es óptimo también ha de ser óptimo el camino de vz al 
vértice v;,, pues si no lo fuese llegaríamos a una contradicción. Si no lo fuese y 
existiera otro camino mejor, incluyendo a éste en el recorrido original 
obtendríamos un camino mejor que el óptimo, lo cual es imposible. Por tanto, se 
cumple el principio de óptimo. 

Planteemos entonces la relación en recurrencia. Para ello, llamaremos D(v,,S) a 
la longitud del camino mínimo que partiendo del vértice v; pasa por todos los 
vértices del conjunto S y vuelve al vértice v;. La solución al problema del viajante 
vendrá dada entonces por D(v¡,V-(v1]): 


D(v,V —(v,)= MiniL(v,,v;) + D(v, Y —(V,V,F 3 
Generalizando para comenzar el recorrido desde cualquier vértice: 
D(v,V)= Min íL(v,,v,) + D(v,,V —v,)j 
¡eV ¡eV 
DO, 0)=L(,,v,) paral<i<n 


Obsérvese la diferencia que existe entre la estrategia de este algoritmo y los que 
tratamos de diseñar siguiendo dos técnicas ávidas (ver el problema 4.4 del capítulo 
anterior). En los algoritmos ávidos se ha de escoger una de las posibles opciones en 
cada paso, y una vez tomada —o descartada—, ya no vuelve a ser considerada nunca. 
Son algoritmos que no guardan “historia”, y por tanto no siempre funcionan. Sin 
embargo, en la Programación Dinámica la solución al problema total se va 
construyendo de otra forma: a partir de las soluciones óptimas para problemas más 
pequeños. 

No obstante, el diseño aquí realizado tiene un serio inconveniente: su 
implementación utilizando una estructura de datos que permita reutilizar los 
cálculos. Tal estructura debería contener las soluciones intermedias necesarias para 
el cómputo de D(v,,V-(v,)), pero estas son demasiadas. 
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En efecto, la tabla debe tener n filas, y 2” columnas, pues éste es el cardinal de 
las partes del conjunto V, que son todas las posibilidades que puede tomar el 
segundo parámetro de D en su definición. 


Por tanto, sí existe una solución al problema del viajante utilizando 
Programación Dinámica, pero no sólo no consigue mejorar la eficiencia de su 
versión clásica mediante Vuelta Atrás (véase el siguiente capítulo), sino que 
tampoco ofrece una mejora en cuanto a la simplicidad de su implementación. 


5.15 HORARIOS DE TRENES 


Una compañía de ferrocarriles sirve n estaciones S,...,S, y trata de mejorar su 
servicio al cliente mediante terminales de información. Dadas una estación origen 
So y una estación destino Sy, un terminal debe ofrecer (inmediatamente) la 
información sobre el horario de los trenes que hacen la conexión entre S, y Sy y 
que minimizan el tiempo de trayecto total. 


Necesitamos implementar un algoritmo que realice esta tarea a partir de la tabla 
con los horarios, suponiendo que las horas de salida de los trenes coinciden con las 
de sus llegadas (es decir, que no hay tiempos de espera) y que, naturalmente, no 
todas las estaciones están conectadas entre sí por líneas directas; así, en muchos 
casos hay que hacer transbordos aunque se supone que tardan tiempo cero en 
efectuarse. 


Solución (O) 


Llamaremos 7(i,j, V) al tiempo del trayecto mínimo para ir de la estación de origen ¡ 
a la estación destino j, pudiendo utilizar como estaciones intermedias las 
contenidas en el conjunto Y, y llamaremos L(i,7) al tiempo del trayecto directo de ¡ 
a j, siendo co si esta conexión no existe. De forma análoga a como razonábamos 
para el problema de los embarcaderos sobre el río (apartado 5.6), podemos plantear 
la solución al problema mediante una ecuación en recurrencia: 


T(G¡,jV)= Min ÁL(, 7),T(k, j,V —k)+LG, 0) 


keV kzi kz j 


Esta ecuación trata de comprobar si es más beneficioso ir de forma directa o a 
través de cada uno de los posibles caminos. Esto hay que hacerlo para cada par (1,7) 
desde 1 hasta n. Obsérvese cómo llegamos a ella por el principio de óptimo pues 
cualquier subtrayecto de un trayecto óptimo a de ser, a su vez, óptimo. 


Podemos representar nuestro problema mediante un grafo siendo las estaciones 
los vértices del grafo y las aristas las conexiones entre dos estaciones, pudiendo no 
existir si no hay trayecto directo entre ambas. La solución al problema se puede 
alcanzar resolviendo el algoritmo de Dijkstra para cada vértice del grafo, puesto 
que tal algoritmo calcula los caminos mínimos desde un único origen hasta los 
demás vértices en un grafo. 


Respecto a su complejidad, ésta coincide con la del algoritmo de Dijkstra, que 
es aceptable para un número n de estaciones razonablemente grande. 


202 TÉCNICAS DE DISEÑO DE ALGORITMOS 


5.16 LA MOCHILA (0,1) 


En el apartado 4.7 del capítulo anterior se planteó el problema de la Mochila (0,1), 
que consistía en decidir de entre n objetos de pesos pj, P>,..., Pn y beneficios 
bi, b»,..., Dn, cuáles hay que incluir en una mochila de capacidad M sin superar 
dicha capacidad y de forma que se maximice la suma de los beneficios de los 
elementos escogidos. Los algoritmos ávidos planteados entonces no conseguían 
resolver el problema. Nos cuestionamos aquí si este problema admite una solución 
mediante Programación Dinámica. 


Solución (8) 


Para encontrar un algoritmo de Programación Dinámica que lo resuelva, primero 
hemos de plantear el problema como una secuencia de decisiones que verifique el 
principio de óptimo. De aquí seremos capaces de deducir una expresión recursiva 
de la solución. Por último habrá que encontrar una estructura de datos adecuada 
que permita la reutilización de los cálculos de la ecuación en recurrencia, 
consiguiendo una complejidad mejor que la del algoritmo puramente recursivo. 


Siendo M la capacidad de la mochila y disponiendo de n elementos, llamaremos 
V(i,p) al valor máximo de la mochila con capacidad p cuando consideramos i 
objetos, con 0<p < My 1<i¡S< mn. La solución viene dada por el valor de V(n,M). 
Denominaremos d, d,, ..., d, a la secuencia de decisiones que conducen a obtener 
V(n,M), donde cada d; podrá tomar uno de los valores 1 ó 0, dependiendo si se 
introduce o no el ¡-ésimo elemento. Podemos tener por tanto dos situaciones 
distintas: 


e Que d, = 1. La subsecuencia de decisiones d;, d», ..., d,, ha de ser también 
óptima para el problema V(n-1,M-p,), ya que si no lo fuera y existiera otra 
subsecuencia e, e>, ..., €,-1 Óptima, la secuencia e, €», ..., €n-1, d, también sería 
Óptima para el problema V(n,M) lo que contradice la hipótesis. 


e Que d, = 0. Entonces la subsecuencia decisiones d, d,, ..., d,., ha de ser también 
óptima para el problema V(n-1,M). 


Podemos aplicar por tanto el principio de óptimo para formular la relación en 
recurrencia. Teniendo en cuenta que en la mochila no puede introducirse una 
fracción del elemento sino que el elemento ¡se introduce o no se introduce, en una 
situación cualquiera V(i,p) tomará el valor mayor entre V(i-1,p), que indica que el 
elemento ¡ no se introduce, y V(i-1,p-p,)+b;, que es el resultado de introducirlo y 
de ahí que la capacidad ha de disminuir en p; y el valor aumentar en b,, y por tanto 
podemos plantear la solución al problema mediante la siguiente ecuación: 


0 sii=0yp>0 
V(,p)= 29 sip<0 
MaxiV(i-1, p)V(i-1,p-—p,)+b,j en otro caso. 


Estos valores se van almacenando en una tabla construida mediante el algoritmo: 
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TYPE TABLA = ARRAY[1..n],[O..M] OF CARDINAL; 
DATOS = RECORD peso,valor:CARDINAL; END; 
TIPOOBJETO = ARRAY [1..n] OF DATOS; 


PROCEDURE Mochila (i,p:CARDINAL; VAR obj: TIPOOBJETO) : CARDINAL; 
VAR elem,cap: CARDINAL; V: TABLA; 
BEGIN 
FOR elem:=1 TO i DO 
Vl[elem,0] :=0; 
FOR cap:=1 TO p DO 
IF (elem=1) AND (cap<obj[1].peso) THEN 
V[elem,cap] :=0 
ELSIF elem=1 THEN 
V[elem,cap] :=0bj[1] .valor 
ELSIF cap<objlelem] .peso THEN 
V[elem,cap] :=V[elem-1,cap] 
ELSE V[elem,cap] := 
Max2(V [elem-1,cap] ,objlelem] .valor+V[elem-1,cap-obj[elem] .peso]) 
END 
END 
END; 
RETURN V[i,p] 
END Mochila; 


El problema se resuelve invocando a la función con ¡=n, p=M. La complejidad 
del algoritmo viene determinada por la construcción de una tabla de dimensiones 
nxM y por tanto su tiempo de ejecución es de orden de complejidad O(1M). La 
función Max2 es la que calcula el máximo de dos valores. 

Si además del valor de la solución óptima se desea conocer los elementos que 
son introducidos, es decir, la composición de la mochila, es necesario añadir al 
algoritmo la construcción de una tabla de valores lógicos que indique para cada 
valor El i,¡] si el elemento ¿ forma parte de la solución para la capacidad ¡ o no: 


TYPE ENTRAONO = ARRAY[1..n],[O..P] OF BOOLEAN; 


PROCEDURE Max2especial(x,y:CARDINAL;VAR esmenorx: BOOLEAN) : CARDINAL; 
BEGIN 
IF x>y THEN 
esmenorx:=FALSE; 
RETURN x 
ELSE 
esmenorx:=TRUE; 
RETURN y 
END 
END Max2especial; 
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PROCEDURE Mochila2(i,p:CARDINAL;obj:TIPOOBJETO;VAR E:ENTRAONO) 
: CARDINAL; 
VAR elem,cap: CARDINAL; V: TABLA; 
BEGIN 
FOR elem:=1 TO i DO 
Vl[elem,0] :=0; 
FOR cap:=1 TO p DO 
IF (elem=1) AND (cap<obj[1].peso) THEN 
V[elem,cap] :=0; 
E[elem,cap] :=FALSE 
ELSIF elem=1 THEN 
V[elem,cap] :=0bj [1] .valor; 
E[elem,cap] : =TRUE 
ELSIF cap<objlelem] .peso THEN 
V[elem,cap]:=V[elem-1,cap]; 
E[elem,cap] :=FALSE 
ELSE V[elem,cap]:=Max2especial (V[elem-1,cap], 
objlelem] .valor+V[elem-1,elem-obj [elem] .peso] ,Elelem,capl); 
END 
END 
END; 
RETURN V[i,p] 
END Mochila2; 


Por otra parte, es necesario construir un algoritmo que interprete los valores de 
esta tabla para componer la solución. Esto se realizará recorriéndola en sentido 
inverso desde los valores ¡ = n, ¡ = M hasta i¡ = 0, ¡ = 0, mediante el siguiente 
algoritmo: 


TYPE SOLUCION = ARRAY[1..n] OF CARDINAL; 


PROCEDURE Componer(VAR sol: SOLUCION); 
VAR elem,cap: CARDINAL; 
BEGIN 
FOR elem:=1 TO n DO (* inicializa solucion *) 
sol [elem] :=0 
END; 
elem:= n; cap:= M; 
WHILE (elem<>0) AND (cap<>0) DO 
IF entralelem,cap] THEN 
sol [elem] :=1; 
cap:=cap-obj[elem] .peso 
END; 
DEC (elem) 
END 
END Componer'; 
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5.17 LA MOCHILA (0,1) CON MÚLTIPLES ELEMENTOS 


Este problema se basa en el de la Mochila (0,1) pero en vez de existir n objetos 
distintos, de lo que disponemos es de » tipos de objetos distintos. Con esto, de un 
objeto cualquiera podemos escoger tantas unidades como deseemos. 


Este problema se puede formular también como una modificación al problema 
de la Mochila (0,1), en donde sustituimos el requerimiento de que x=06 x= 1, 
por el que x; sean números naturales. Como en el problema original, deseamos 
maximizar la suma de los beneficios de los elementos introducidos, sujeta a la 
restricción de que éstos no superen la capacidad de la mochila. 


Solución (O) 


Para encontrar un algoritmo de Programación Dinámica que resuelva el problema, 
primero hemos de plantearlo como una secuencia de decisiones que verifiquen el 
principio del óptimo. De aquí seremos capaces de deducir una expresión recursiva 
de la solución. Por último habrá que encontrar una estructura de datos adecuada 
que permita la reutilización de los cálculos de la ecuación en recurrencia, 
consiguiendo una complejidad mejor que la del algoritmo puramente recursivo. 


Con esto en mente, llamaremos V(i,p) al valor máximo de una mochila de 
capacidad p y con i tipos de objetos. Iremos decidiendo en cada paso si 
introducimos o no un objeto de tipo ¿. Por consiguiente, para calcular V(i,p) existen 
dos opciones en cada paso: 


e No introducir ninguna unidad del tipo ¿, con lo cual el valor de la mochila V(i,p) 
es el calculado para V(i-1,p). 


e Introducir una unidad más del objeto ¡ lo cual indica que el valor de V(1,p) será 
el resultado obtenido para V(i,p-p;) más el valor del objeto v;, con lo cual se 
verifica que V(i,p) = V(i,p-p;) + b;. 


Esto nos permite establecer la siguiente relación en recurrencia para V(i,p): 


0 si p=0 
V(i,p)= (p=p,)b, si ¿=1 
V(i=1, p) si p<p, 


MaxV(i¡—1L p)V(i,p—p,)+b,) en otro caso. 


Utilizaremos una matriz nxM para almacenar los valores de Y que vayamos 
obteniendo y así no repetir cálculos. El algoritmo que resuelve el problema es el 
siguiente: 


TYPE TABLA = ARRAY [1..n],[O0..M] OF CARDINAL; 
DATOS = RECORD peso,valor:CARDINAL; END; 
TIPOOBJETO = ARRAY [1..n] OF DATOS; 
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PROCEDURE Mochila3(i,p:CARDINAL;VAR obj:TIPOOBJETO) : CARDINAL; 
VAR elem,cap: CARDINAL; V: TABLA; 


BEGIN 
FOR elem:=1 TO i DO (* condiciones iniciales *) 
Vl[elem,0] :=0 
END; 


FOR cap:=1 TO p DO 
V[1,capl:=(cap DIV obj[1].peso)*obj [1] .valor 
END; 
FOR elem:=2 TO num DO 
FOR cap:=1 TO p DO 
IF (cap < plelem]) THEN V[elem,cap]:=V[elem-1,capl] 
ELSE V[elem,cap]:= Max2(V[elem-1,cap], 
(V[elem,cap-obj [elem] .peso]+0bj [elem] .valor)) 
END 
END 
END; 
RETURN V[i,p] 
END Mochila3; 


La complejidad del algoritmo es la que corresponde a la construcción de la 
tabla, es decir O(1M). 


5.18 LA MULTIPLICACIÓN ÓPTIMA DE MATRICES 


Necesitamos calcular la matriz producto M de n matrices dadas M = M¡M»...M,, 
minimizando el número total de multiplicaciones escalares a realizar. Este 
problema ya fue planteado en el capítulo anterior, en donde vimos que los 
algoritmos ávidos presentados no encontraban solución en todos los casos. 

Nos preguntamos ahora si existe un algoritmo que lo resuelva utilizando la 
Programación Dinámica. 


Solución (O) 


En primer lugar, vamos a suponer como hacíamos en el capítulo anterior que cada 
M, es de dimensión d; ¡xd; (1 < i < n), y por tanto realizar la multiplicación M; M;, 
va a requerir un total de d; ¡d,d;+, operaciones. Llamaremos M(i,j) al número 
mínimo de multiplicaciones escalares necesarias para el cómputo del producto de 
M; M;+1 ... M,con 1 < ¡<< n. Por consiguiente, la solución al problema planteado 
coincidirá con M(1,n). 

Para plantear la ecuación en recurrencia que define la solución, supongamos que 
asociamos las matrices de la siguiente manera: 


(M; Miz1 ... Mi) (Mir Mio ... Mj). 
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Aplicando el principio de óptimo, el valor de M(i,7) será la suma del número de 
multiplicaciones escalares necesarias para calcular el producto de las matrices 
M; M;1...Mj, que corresponde a M(i,k), más el número de multiplicaciones 
escalares para el producto de las matrices My+1My+> ...M;, que es M(k+1,7), más el 
producto que corresponde a la última multiplicación entre las matrices de 
dimensiones 
(d; 1 di) y (d; dj) es decir, d; :did; En consecuencia, el valor de M(i,j) para la 
asociación anteriormente expuesta viene dado por la expresión: 


M(ij) = Mk) + MU] + dnde, 


Pero k puede tomar cualquier valor entre i y j-1, y por tanto M(1,7) deberá 
escoger el más favorable de entre todos ellos, es decir: 


MG, J)=MiniMGk)+ MU+1,J)+d, dd, 
¡<k<j 


lo que nos lleva a la siguiente relación en recurrencia: 
0 sii=j 
M(i,j)= MiníM(i,k) + M(k+1,j)+ d, dd y en otro Caso. [5.4] 
¡<k<j 


Para resolver tal ecuación en un tiempo de complejidad polinómico es necesario 
crear una tabla en la que se vayan almacenando los valores M(1,) (1 <i<j¡<mnm) y 
que permita a partir de las condiciones iniciales reutilizar los valores calculados en 
los pasos anteriores. 


Esta tabla se irá rellenando por diagonales sabiendo que los elementos de la 
diagonal principal son todos cero. Cada elemento M[i,j] será el valor mínimo de 
entre todos los pares (M[1,k] + M[K+1,j]) señalados con la línea de doble flecha en 
la siguiente figura, más la aportación correspondiente a la última multiplicación 
(d; ¡dd;). Los valores que requiere el cálculo de M[i,j] y que el algoritmo reutiliza 
para conseguir una tiempo de ejecución aceptable se encuentran sombreados: 


Al ir rellenando la tabla por diagonales se asegura que esta información 
(M[¡,k],M]x+1,7]) está disponible cuando se necesita, pues cada M]i,j] utiliza para 
su cálculo todos los elementos anteriores de su fila y todos los de su columna por 
debajo suya. 
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Rellenada la tabla, la solución la podemos encontrar en el extremo superior 
derecho, que nos indica el número de multiplicaciones escalares buscado, M] 1,n]. 

Si además queremos conocer cómo es la asociación que corresponde a este 
óptimo es necesario conservar para cada elemento M]i,¡] el valor de k para el cual 
la expresión M[i,k] + M[k+1,5] + d; dd, es mínima, construyendo otra tabla que 
denominamos Factor. 

El siguiente algoritmo Matriz es de creación de las tablas M y Factor. En la 
tabla M se almacenan los valores del número mínimo de multiplicaciones y en la 
tabla Factor la información necesaria para construir la asociación óptima. 


TYPE MATRIZ = ARRAY [1..n],[1..n] OF CARDINAL; 
ORDEN = ARRAY [O..n] OF CARDINAL; (* dimensiones *) 


PROCEDURE Matriz(VAR d:ORDEN;n:CARDINAL;VAR M,Factor:MATRIZ) ; 
VAR i,diagonal: CARDINAL; 


BEGIN 
FOR i:=1 TO n DO 
M[i,i]:=0 
END; 


FOR diagonal:=1 TO n-1 DO 
FOR i:=1 TO n-diagonal DO 
M[i,i+diagonal]:= 
Minimo(d,M,i,i+diagonal,Factor[i,i+diagonal]); 
END 
END 
END Matriz; 


La función Minimo es la que calcula el mínimo de la expresión en recurrencia 
[5.4], y devuelve no sólo el valor de este mínimo, sino el valor de k para el que se 
alcanza (mediante el parámetro k1): 


PROCEDURE Minimo(VAR d:ORDEN;VAR M:MATRIZ;i,j:CARDINAL; 
VAR k1:CARDINAL) : CARDINAL; 
VAR aux,k,min: CARDINAL; 
BEGIN 
min:=MAX (CARDINAL) ; 
FOR k:=i TO j-1 DO 
aux:=M[i,k]+M[x+1,31+d[i-1]*d[x1x*d[3]; 
IF aux<min THEN min:=aux; k1:=k END 
END; 
RETURN min 
END Minimo; 
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Observando el procedimiento Matriz vemos que existe un bucle externo que se 
repite desde diagonal = 1 hasta n — 1, y en su interior un bucle dependiente de la 
iteración estudiada y del valor de diagonal, y que se ejecuta desde 1 hasta 
n — diagonal. En el interior de este bucle hay una llamada al procedimiento Minimo 
que tiene una complejidad del orden del valor de la diagonal, y por tanto el tiempo 
de ejecución del algoritmo es: 


n-—1 n-1 n—1 
A (n— diagonal)diagonal = n > diagonal — > diagonal? = (n* =n)/6 


diagonal=1 diagonal=1 diagonal=1 


por lo que concluimos que su complejidad temporal es de orden O(n?). Por otro 
lado, la complejidad espacial del algoritmo es de orden O(m5). 

En caso que deseemos reconstruir la solución a partir de la tabla Factor, el 
siguiente procedimiento muestra por pantalla la forma de multiplicar las matrices 
para obtener ese valor mínimo: 


PROCEDURE EscribeOrden(VAR Factor:MATRIZ;i,j:CARDINAL); 
VAR k:CARDINAL; 
BEGIN 
IF i=j THEN 
WrStr('M?); 
WrCard(i,0) 
ELSE 
k:=Factor[i,j]; 
WrStr(“(); 
Escribe0rden(Factor,i,k); 
WrStr('x>); 
Escribe0rden(Factor,k+1,j3); 
WrStr()?) 
END 
END EscribeUOrden; 


El algoritmo aquí presentado es capaz de encontrar el valor de la solución 
óptima junto con una de las formas de obtenerla. Sin embargo, existen casos en 
donde puede haber más de una forma de multiplicar las matrices para obtener el 
valor óptimo, como muestra el siguiente ejemplo. 


Sean las matrices M¡ (10x10), M> (10x50) y M3 (50x50). Existen dos formas de 
asociarlas para multiplicarlas, y en ambos casos obtenemos: 


(M,M>,)M5 = 10:10:50 +10-50-50 = 30000 
M¡(M,M5 ) = 10-50-50 +10-10-50 = 30000 
Es posible modificar el algoritmo para que encuentre todas las soluciones que 


llevan al valor óptimo, y para esto es suficiente valerse de la matriz Factor, si bien 
esta modificación reviste poco interés desde un punto de vista práctico. 
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